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Collective Kondo effect in the Anderson-Hubbard lattice 
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Abstract. The periodic Anderson model is extended by switching on a 
Hubbard U for the conduction electrons. We use the Gutzwiller variational 
method to study the nearly integral valent limit. The lattice Kondo energy 
contains the ^/-dependent chemical potential of the Hubbard subsystem in 
the exponent, and the correlation-induced band narrowing in the prefactor. 
Both effects tend to suppress the Kondo scale, which can be understood to 
result from the blocking of hybridization. At half-filling, we find a Brinkman- 
Rice-type transition from a Kondo insulator to a Mott insulator. 
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Recently, there has been a great deal of interest in many-band models 
which describe hybridization between different subsystems of strongly cor- 
related electrons. Schork and Fulde introduced the problem of the Kondo- 
compensation of a magnetic impurity by a correlated band Q. The cor- 
responding periodic model can be envisaged either as an array of localized 
spins coupled to a Hubbard band (the Kondo-Hubbard lattice ||), or as 
a periodic system of Anderson impurities hybridized with a Hubbard band 
(the Anderson-Hubbard lattice |J). A particular version of the latter model 
was suggested for the remarkable correlated semiconductor FeSi |Q . We note 
that a closely related class of two-band models is the subject of intense inves- 
tigation in connection with the resurgent interest in metallic ferromagnetism 

ill- 

In contrast to their single-impurity counterparts, the Anderson-Hubbard 
(AH) and Kondo-Hubbard (KH) lattice models are not equivalent 0, [?]]. 
In fact, they tend to lead to different predictions as to the effect of the 
conduction band Hubbard U on the Kondo energy A strong enhancement 
of the Kondo scale is found for the KH model, both for the single impurity 
|| and the lattice ||. As for AH models, there is at least one contribution 
which indicates the reduction of the Kondo scale with increasing U for the 
impurity [l], ||, and there is an argument for a reduced scale in the lattice 
[[|. The latter finding may be understood as the consequence of the blocking 
of hybridization processes. 
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We consider the periodic Anderson-Hubbard model which describes an 
array of strongly correlated /-sites hybridized with a moderately strongly 
interacting c?-band: 

k,CT j,<7 j 

+^E"/t"S - «E(/jU» + <*U0 (i) 

where h$ a = fj a fj a , etc., the k are wave vectors, and the j are site indices. 
The <i-bandwidth is W. In what follows, we take the strongly asymmetric 
Anderson model with Uf — *■ oo and the /-level e/ < sufficiently deep- 
lying so that we are in the Kondo limit: 1 — rif <C 1 where the /-valence is 
defined as rif = (^2 a n^ a ). The total electron density (per site, for one spin) 
is n = (X)j,<T^j(7 + Sj,o-"jcr)/2-k, where L is the number of lattice sites. We 
will assume 1/2 < n < 1, so that there are enough electrons to fill at least 
the /-levels, and the <i-band filling is variable up to half-filling. 

We use the Gutzwiller variational method, generalizing a previous treat- 



ment of the periodic Anderson model ||10[ . We postulate a Gutzwiller- 
projected hybridized band ground state 

i*> = % • % ■ nriK/L + mljio) (2) 

k cr 

where the mixing amplitudes u^/v^ are treated as independent variational 
parameters. The Gutzwiller projector for the G?-electrons is 

Pg = W-(1-v)KA] ( 3 ) 



where the variational parameter t] is controlled by Ud- For the /-electrons the 
full Gutzwiller projection is taken Pq = Yl s [l—h^h^]. Here we consider only 
non- magnetic solutions corresponding to a mass-enhanced metal at n < 1, 
or a renormalized- hybridization-gap insulator at n — 1. 

Omitting the details of the derivation (which relies on the Gutzwiller 
approximation), we quote our essential results. The expression for the opti- 
mized total energy density is 

£ = e f -n d (l-nO d )g d + Uy d 

The first line gives the energies of decoupled /-, and ci-electrons. The coupling 
between the two subsystems is described by the last term which we identify 
as the Kondo energy of the Anderson-Hubbard model. The same Kondo 
scale can be identified in the deviation of the valence from 1 

l - nf = WJW e ^\ WW) r () 

In the above equations, n d = n — 1/2 is the v = value of the conduction 
electron density (per spin). v\ is double occupation, is the q- factor (the 
renormalization factor of the hopping amplitude) taken with n d , and [io(Ud) 
is the chemical potential of the Hubbard band calculated (in the Gutzwiller 
approximation) for band filling n d . 

Eqs. (|) and (J4|) yet formal results: they express the solution in 

terms of the optimized parameters of the Hubbard subsystem. There is no 



closed-from result away from half-filling; the results of sma\\-Ud, and of large- 
Ud, expansions are given in |||. Here we discuss first the general features, 
and then the special case of half-filling. 

The conduction band Ud appears in the characteristic Kondo scale in two 
ways: through the prefactor and the exponent. The prefactor describes the 
correlation-induced narrowing of the ci-band. The exponential factor can be 
written as exp (— W/J c s) where we introduced the effective Kondo coupling 

Av 2 

Jes = fTT v • (6) 

It differs from the Anderson lattice result |K| [nj inasmuch as /xo(0) = 
(n d — 1/2) W is replaced by the {/^-dependent chemical potential Ho{Ud)- 
Comparing to the single impurity result Q, we call attention to the extra 
factor of 2 in the numerator: the 'lattice enhancement of the Kondo effect' 
[ PH works also in the present model. 

Due to the combined effect of the prefactor and the exponent, the lattice 
Kondo energy is a decreasing function of Ud- To make this finding plausible, 
one can argue that increasing Ud suppresses charge fluctuations, and therefore 
blocks the hybridization processes which give rise to the Kondo coupling, 
and thus to the Kondo effect. The blocking of hybridization is carried to an 
extreme in the case of exact half-filling, as we discuss below: 

According to a usual Luttinger's theorem argument, the half-filled case 
n = 1 belongs to a non-magnetic insulator. The well-known qd = 1 — 
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(Ud/ZW) 2 holds right up to Ud = 2W where the <i-band undergoes a Brinkman- 



Rice transition |L2|]. However, in contrast to the usual one-band Brinkman- 
Rice transition which is a metal-insulator transition, here we have to do with 
an insulator-insulator transition. For Ud < 2W, the system is a new kind of 
Kondo insulator with a renormalized hybridization gap Q. For Ud > 2W, the 
system is integral valent (jif = 1), the Kondo effect is completely quenched, 
and we should speak about a Mott insulator. This is certainly what we ex- 
pect within the Gutzwiller approximation. Though the well-known criticism 
concerning the complete suppression of polarity fluctuations can be raised 
here, we note that up to the Brinkman-Rice transition, the Gutzwiller re- 
sult is in a surprisingly good agreement with the predicted high-dimensional 
behaviour of the Hubbard model The variationally derived lattice en- 
hancement of the Kondo scale is at least qualitatively confirmed by a similar 
calculation for the Anderson lattice [|14j]. Therefore we expect that our re- 



sults which combine Hubbard and Anderson lattice features, have a similar 
justification for Ud < 2W. For Ud > 2W, a treatment based on replacing 
the Hubbard band with a t-J |T5j or a Heisenberg model [ITBJ may be more 



appropriate. 

Our trial state is the lattice generalization of a lowest-order Varma-Yafet 



state |L7| which is a potential source of shortcomings. Work on the impurity 
case is in a more advanced state and indicates that processes which enhance 
the Kondo scale, may predominate || ||. It has been argued || that to 
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capture this effect variationally, one should consider states which are higher 
in the Varma-Yafet hierarchy. If a similar statement proves true for the 
lattice, our work served to demonstrate the limitations of the Ansatz (|2|). 
However, we believe that the lattice case may be different, for the following 



reason: By postulating a large Fermi surface |18|], states from the entire 
Brillouin zone are used, and therefore (§) does contain contributions from 
electron-hole pairs with respect to the small (conduction electron) Fermi 
surface. 

Acknowledgement. P. F. gratefully acknowledges financial support by 
the Hungarian National Science Research Foundation grant OTKA T-014201. 

References 

[1] T. Schork and P. Fulde, Phys. Rev. B 50 (1994) 1345. 

[2] N. Shibata et al., Phys. Rev. B 53 (1996) R8828. 

[3] K. Itai and P. Fazekas, Phys. Rev. B 54 (1996) R752. 

[4] C. Fu and S. Doniach, Phys. Rev. B 51 (1995) 17439. 

[5] Hal Tasaki, Phys. Rev. Lett. 75 (1995) 4678. 

[6] K. Penc, H. Shiba, F. Mila, and T. Tsukagoshi: Phys. Rev. B 54 (1996) 
4056. 



7 



[7] Y.M. Li, Phys. Rev. B 52 (1995) 6979. 

[8] G. Khaliullin and P. Fulde, Phys. Rev. B 52 (1995) 9514. 

[9] T. Schork: Phys. Rev. B 53 (1996) 5626. 

[10] P. Fazekas and B.H. Brandow, Phys. Scr. 36 (1987) 809. 

[11] T.M. Rice and K. Ueda, Phys. Rev. B 36 (1986) 6420. 

[12] W.F. Brinkman and T.M. Rice, Phys. Rev. B 2 (1970) 4302. 

[13] G. Moeller et al, Phys. Rev. Lett. 74 (1995) 2082. 

[14] M. Jarrell, Phys. Rev. B 51 (1995) 7429. 

[15] S. Moukouri et al., Phys. Rev. B 53 (1996) R488. 

[16] J. Igarashi et al., Phys. Rev. B 51 (1995) 5814; ibid 52 (1995) 15966. 

[17] CM. Varma and Y. Yafet, Phys. Rev. B 13 (1976) 2950. 

[18] H. Shiba and P. Fazekas, Progr. Theor. Phys. Suppl. 101, 403 (1990). 



8 



